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Abstract 

m 

An asymptotic technique is presented to characterize the bits/symbol achievable on a representative wireless 
I link in a spatially distributed network with active interferers at correlated positions, N receive diversity branches, 

■ and linear Minimum-Mean-Square-Error (MMSE) receivers. This framework is then applied to systems including 
I analogs to Matern type I and type II networks which are useful to model systems with Medium-Access Control 

■ (MAC), cellular uplinks with orthogonal transmissions and frequency reuse, and Boolean cluster networks. It is 
I found that for our network models, with moderately large N, the correlation between interferer positions does 

■ not significantly influence the bits/symbol resulting in simple approximations for the data rates achievable in such 



(N 



networks which are known to be difficult to analyze and for which few analytical results are available. 

Index Terms 



^ • Matern, MIMO, MMSE, Cellular, Reuse. 

I. Introduction 

^ ' In systems with multiple diversity branches such as multi-antenna systems, the diversity can be used 
for interference mitigation which can improve spatial reuse, and increase aggregate data rates in networks. 
O Since signal and interference strengths are influenced by relative node positions, the spatial distribution 
^ of nodes strongly influences the performance of such systems. Most works on wireless networks with 
^ spatially distributed nodes assume that the locations of active nodes are independent, a good model for 
networks with uncoordinated transmissions (e.g. [H], [IS, ill, Q). In many situations however, active 
^ nodes may be spatially correlated, e.g. in networks with coordinated transmissions such as networks with 
>• Medium-Access-Control (MAC), systems with significant physical constraints on how close nodes may 
^ get to each other such as vehicular networks, networks which employ location-based spatial reuse such 
5-j as cellular networks, and networks with spatially clustered users. The analysis of networks with spatially 
correlated transmitters is well known to be difficullQ. 

In this paper, we develop an asymptotic technique to characterize the Signal-to-Interference-Ratio 
(SIR) of a representative link with multiple receive-diversity branches employing a linear Minimum- 
Mean-Square-Error (MMSE) receiver, which is the linear receiver that maximizes the SIR if noise is 
negligible compared to interference. The SIR is a useful metric in dense wireless networks which are 
interference-limited and directly influences achievable data rates if Gaussian transmit signals are assumed^. 
The representative link is assumed to operate in a network with co-channel transmitters at correlated spatial 
positions. The asymptotic regime we consider is when the number of diversity branches N and potential 
transmitting nodes n (to be defined in Section II) go to infinity with a fixed positive ratio. We find that 

This research was supported in part by the National Science Foundation under Grant CCF-1 117218. Portions of this work were presented 
at the 2012 International Symposium on Information Theory (ISIT). 

'We use the term spatially correlated in this work to refer to networks where the position in space of randomly distributed transmitters is 
coiTelated, as compared to systems where channel fading between antennas are correlated 

^Without the assumption of Gaussian transmit signals, we can apply a correction factor to the SIR to estimate the bits/symbol (6). 
Moreover, since Gaussian interference maximizes the interference entropy, the bits/symbol computed by assuming Gaussian transmissions 
can be interpreted as a lower bound on the bits/symbol with non-Gaussian transmit signals. 
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an appropriately normalized version of the SIR, and under the assumption of Gaussian transmissions, 
the achievable bits/symbol converge in probability to deterministic quantities if the correlation between 
active transmitters satisfies an asymptotic independence property. The normalized mean bits/symbol is also 
found to converge to the same quantity. The asymptotic mean bits/symbol for different spatial correlation 
models that satisfy the asymptotic independence property, which includes independent transmissions, 
take a similar form, but with different effective densities of nodes. Thus, the asymptotic bits/symbol 
for systems with "large" numbers of diversity branches and correlated transmissions (satisfying the 
asymptotic independence property), is approximated well by the asymptotic bits/symbol for a network 
with independent transmissions but equivalent density of active nodes. This framework is applied to two 
hard-core models which are close analogs to the Matem type-I and type-II networks commonly used 
to model Carrier-Sense-Multiple-Access (CSMA) networks, cellular networks with frequency reuse, and 
a Boolean cluster model whereby nodes are active if they are within a certain distance of randomly 
distributed cluster center, which illustrates the applicability of this framework to a wide range of systems. 

The majority of works which consider spatial correlation models which admit hard-core networks deal 
with systems that have no diversity (e.g. single antenna/node) and use asymptotic techniques or other 
approximations as hard-core networks are known to be particularly difficult to analyze [7|. For instance, 
Haenggi, in [7J, analyzes Matem Type-I and II hard-core networks with single- antenna nodes by finding 
bounds on the mean interference. Ganti, Andrews and Haenggi in [8J, and Giacomelli, Ganti and Haenggi 
in im, analyze networks with single- antenna nodes at correlated positions (including hard-core networks) 
using asymptotic techniques. In particular they consider the high SIR and reliability regimes respectively, 
where the density of active nodes is taken to zero. The few prior works that consider multi-antenna systems 
in hard-core networks (e.g. [lOJ) approximate the positions of the active nodes by a non-homogenous 
Poisson process by reducing the density of active nodes in the vicinity of a test receiver, but assuming 
that node positions are uncorrected. In contrast to these works, the hard-core models we analyze (which 
are close analogs to the Matem type-I and II systems) explicitly consider spatially correlated transmitting 
nodes, with receive diversity. The complexity caused by the spatial correlation of active nodes is handled 
by the asymptotic analysis. 

The analysis of cellular networks with regular cells, even for the simpler case of universal frequency 
reuse, is known to be very difficult (as noted in ["61 and flT]), and are usually analyzed by Monte-carlo 
simulation such as in [,12,1 and other references given in [1 1] . In a recent work [13] , we consider the uplink 
of a cellular system in a similar asymptotic regime but with the active nodes at independent locations 
on the plane. In contrast, the framework developed in this work is applicable to cellular networks which 
use orthogonal multiple-access, whereby at most one user is allowed to transmit in a given time slot per 
cell, and nearby-cells occupy different frequency bands according to some frequency reuse pattem, which 
results in spatially correlated transmitters. 

Of the network models with correlated active users, clustered wireless networks have proven to be 
somewhat tractable compared to the hard-core networks and frequency-reuse systems discussed above. In 
particular, networks with Poisson clusters, whereby users are distributed independently in clusters, with 
the cluster centers forming a homogenous Poisson Point Process (PPP) on the plane have been analyzed 
in works such as [14] and [|T5]| . The main feature of Poisson cluster processes that makes them amenable 
to analysis is the high degree of independence inherent in the model, since the locations of users within 
each cluster are independent. Furthermore, the cluster locations are independent as well. The framework 
we develop in this work however, is applicable to a different cluster model, namely the Boolean cluster 
model whereby nodes are distributed independently on the plane, and nodes are actively transmitting if 
they are within a certain distance of at least one cluster center. The cluster centers are distributed randomly 
and independently on the plane. A major difference between the Boolean cluster model and the Poisson 
cluster model is that in the former, the density of active nodes on the plane is either zero or a positive 
constant. In the Poisson cluster model however, the density of nodes on the plane could vary significantly 
as regions where K > 1 clusters overlap say, will have K times the density of nodes compared to regions 
of the plane where there is a single cluster. 
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The asymptotic analyses carried out in this work are vaUdated by Monte Carlo simulations which 

show a close agreement between the asymptotic predictions (particularly the mean bits/symbol) and 
simulated values, even when the diversity order per receiver is only moderately large (e.g. 6). This 
finding is consistent with what has been frequently observed in the random matrix literature, whereby 
the convergence of various quantities to their asymptotic values are rapid. Note that while it would be 
ideal to have an exact statistical characterization of the Signal-to-Interference-plus-Noise Ratio (SINR) 
for finite systems with the forms of spatial correlations that fall under our framework, we note that such 
systems have proven extremely difficult to analyze exactly, even for single antenna systems, as noted in 
the references summarized above and evidenced by the approximations that are frequently used in the 
literature to analyze such networks. Thus, the asymptotic technique developed here can be used to gain 
insight into the performance of networks which are known to be very difficult to analyze, and for which 
few analytical results of any form exist in the literature. 

II. System Model 

Consider a circular network of radius R centered at the origin, with a representative receiver at the 
origin. The representative receiver is in a link with a representative transmitter located at a deterministic 
point Xt which is at distance = \X-r\ from the origin. The representative receiver has N diversity 
branches (e.g. antennas or frequency diversity), with independent flat fading on each diversity branch. In 
addition to the representative transmitter, there are n potential interferers in the network, which if active, 
are co-channel interferers to the representative link. We use the word potential here since some of the 
potential interferers will not be actively transmitting as described later in this section. Let the potential 
interferer locations be denoted Xi, X2, • • • , X„ with the distance of the i-th potential interferer from the 
origin denoted by = The n potential interferers are i.i.d. with uniform probability in the circular 
network such that the area density of potential interferers pp satisfies 

n = n Pp R^. (1) 

The average power (averaged over the fading realizations) received at each diversity branch of the repre- 
sentative receiver from an transmitter at distance fransmitting with power is pi = Pir^°', with a > 2. 
The transmit power of the i-th potential interferer, Pj equals one or zero, and the transmit power of the 
representative transmitter is fixed at unity. Therefore, the vector of transmit powers p = (Pi, P2, ■ ■ ■ , P„) 
controls which of the potential interferers are active, and the representative transmitter is assumed to 
always be active. 

The vector y defined as follows contains the sampled signals at a given sampling time at the N diversity 
branches of the representative receiver. 

n 

-- \ — " -- / 

y = ' StXt + X^r^ ^ giVPiXi + w , (2) 

i=l 

where Xi (xt), is the transmitted symbol of potential interferer i (representative transmitter), and gj G 
^Nxi g C^^^) contains i.i.d., zero-mean, unit variance random variables drawn from a continuous 
distribution, representing the fading between the i-th potential interferer (representative transmitter) and 
the N diversity branches of the representative receiver, w G C^^^ represents independent zero-mean, 
circularly-symmetric Gaussian noise of variance cr^ per complex dimension. In this paper, we shall focus on 
the interference-limited regime and neglect the noise by setting = 0. The active nodes in the network are 
assigned non-zero transmit powers by a function g{-) such that p = g {Xi, • • ■ , X^, Mi, M2, • • ■ , M,n; Xt). 
Ml, ■ ■ ■ , Mrn are auxiliary random variables that are used to select the active transmitters and are defined 
differently according to the specific spatial correlation model. For instance, in the HC-II model these 
variables are used to select which node is active when there are two or more nodes within each others 
hard core, and in the Boolean cluster model, these variables represent the cluster centers. The function g{-), 
thus controls which of the n potential interferers are transmitting, and since it is parameterized by Xt, 
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the transmit powers of the interferers are a function of the locations of all potential transmitters (potential 
interferers and the representative transmitter) in the network. We denote the set of active transmitters as 

T:= {i:P, = l}. 

The asymptotic regime we shall consider is the limit as N, n and R ^ oo, such that n/N = c > 
and pp are constants, and ([T]) holds. For the rest of this paper, whenever any one of the quantities n, N 
or i? — 7- oo, it is assumed that the other two quantities go to infinity as well, such that n/N = c and 
([U) hold. The main results are given in terms of limiting values of a normalized version of the SIR, 
= iV-"/V^SIR, at the output of the MMSE receiver. This normalization is used because the (un- 
normalized) SIR grows with the number of diversity branches as A^"/^. Thus, using the standard formula 
for the SIR at the output of an MMSE receiver, we have 

Pn = N-^ gt. P, rr- g, gt j g^ . (3) 

Define pin = N'^/^P.r^'^, and pj^ = N'^'^P^rJ'' for convenience. The spatial correlation between the 
transmitting nodes needs to satisfy the following asymptotic independence property, 

^ n n 

lim ^ [Pr(pin < x,pjn <x)- Pr(j9i„ < x) Pr(pj„ <x)] = Q. (4) 

1=1 j=l 

If there is symmetry between the indices of the potential interferers, i.e. Pr(pi„ < x) = Pr(pj„ < x) 
and for all i ^ j and k ^ i, Fi^Pm < x,pjn < x) = Fi{pkn < x,pen < x), then dH) directly simplifies to 
the following for i ^ j 

lim Piipin < x,pjn <x)= lim Pr(pi„ < x) Pr(]9j„ < x) . (5) 

Note that Pin is the average received power (if averaged over the fading distribution) from node i scaled 
by A^"/^, as observed at the origin. Furthermore, note that although we refer to dD and dD as asymptotic 
independence, it is a slightly weaker requirement than asymptotic independence as the joint CDFs and 
the product of the marginal CDFs of pin and pjn that appear in dU and dS) are evaluated at a common 
argument x. The asymptotic independence property is of course trivially satisfied if the locations of the 
active nodes are independent. Additionally, define u as the average probability of a potential interferer 
being active in the limit, given by 



1/ = lim 



1 " 

-VPr(Pi = l). (6) 



n— >-oo Tl 

i=l 

We require that u/N > 1, in other words, there are more active interferers than the diversity order at the 
receiver in the limit. 

III. Main Result 

The empirical distribution function (e.d.f.) of a set of random variables, which is used in the main 
results, is defined as the proportion of those random variables that are less than or equal to the argument 
of the e.d.f. For instance, denoting the e.d.f. of pin by Hn{x), we have 



n 

HJx) = -S^l 



1 " 

} = -y]am(a;), (7) 



1=1 i=l 

where is the indicator function, and ain{x) = l{pi„<x}- We can now state the main result. 

Theorem 1: If the asymptotic indepence property defined in dl) holds, as n,N,R — )• oo such that 
n/N = c > 0, and ([U) hold, /^at — )■ /3 in probability where (3 is the real, positive solution to 

Jo 1 + r/3 
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where H{x) is defined as follows 

lim Hn{x) — > H{x) , in probability . (9) 

n—>oo 

The convergence in probability above is guaranteed by (H). Additionally, if the following holds, 

lim Pr(rr°A^f < x\Pi = 1) = lim Pr(r-"A^t < x) , (10) 

n— >oo n— >oo 

then dH]) evaluates to 

?f£ii^J?l)^l^ ^^-Pi'-'P ,,f.fl_i.l_l;2-^;^V (11) 

tt V"/ (« - 2)(c + 7rp/3)^"s: \ a a a 7rpfj + cJ 

where p = ppU, and 2-^i(-, •;•;•) is the Gauss hypergeometric function (e.g. see [16J). 
Proof: Given in Appendix El 

Note that p is the effective density of active transmissions in the limit. Additionally, note that when the 
is large, the SIR is not strongly influenced by the dependence of the transmit powers as /3 is only dependent 
on the marginal probability that a node is selected to transmit. This leads to simple approximations to the 
bits/symbol given later in this section as the complicated dependence between transmit powers becomes 
negligible with large N. 

If we assume that the locations of the potential interferers are independent, then the asymptotic 
independence property is trivially satisfied, which means Theorem 1 holds. Hence, the normalized SIR 
in systems with correlated active transmitters, (but satisfying the asymptotic independence property), 
converges to the same limit as a network with independent transmitters with equal density of transmitting 
nodes. 

Assuming that all nodes use Gaussian code books and recalling that the thermal noise is negligble, we 
estimate the bits/symbol using the Shannon formula, Cn = log2(l + SIR), where the subscript is used 
to emphasize the dependence of the bits/symbol on the number of receiver degrees of freedom. By the 
continuity of the log function, 

|log2(l + SIR) - log2 (1 + r^"Art/3) | ^ o 

in probability as A^ — )• oo in the manner of Theorem 1 . Furthermore, applying the bounded convergence 
theorem in the manner of UJJ, we can show that 

E [log2(l + SIR)] - log2 (1 + r^^ATf /3) ^ 0. 

Hence, log2 (l + r^'^N^ (3) is a good approximation for the mean bits/symbol when the diversity order 
A^ is large. 

In general, the solution for /3 in ([8]) or (fTTT) has to be found numerically. However, if (fTOl ) holds, and the 
diversity order at the receiver A^ is much smaller than the number of potential interferers in the network 
n, c is large, and we can find a simple approximation for (3 that exposes the dependency of the SIR on 
system parameters. Note that the hypergeometric function in the second term on the RHS of (fTTI) has 
positive parameters implying that it is a power series with positive coefficients. Hence, the second term 
on the RHS of (fTTI) is non-negative and is an increasing function of its argument, and thus is bounded 
from above by 

^'^"'"^^ ,,f,fl_i,l_l;2-l;l)=J^M!i4cscfH^) (12) 



a 



{a-2){c + 7rp(3y-^ \ a a a J a{c + Tip(3Y 

Since a; > 2, as c — )• oo, the previous expression goes to zero. Thus the second term on the RHS of (fTTI) 
goes to zero as well. Thus, if the number of potential interferers in the network n, greatly exceeds the 
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diversity order at the receiver, i.e. c is very large, the second term on the RHS of (fTTT ) is small and {3 can 
be approximated as 



/3 





a 1 

















(13) 



Thus, when and n are both large but n N, the normalized SIR (3^ can be approximated by 
p. Removing the normalization by A^"/^ and applying the Shannon formula leads to the following 
approximation for the bits/symbol and its mean 



C ^ E[C] ^ Iog2 1 + 



^ o 2 2 sm — 



(14) 



IV. Application to Matern-like Hard Core Models 

We apply the framework developed in Section III to two hard-core models, HC-I and HC-II which 
are analogous to the Matem Type I and II models commonly used to model hard-core processes [18|. 
Hard-core models are simple models for active nodes in systems with MAC protocols where nodes that 
are close to each other are not allowed to transmit simultaneously. For instance, the Matem Type II model 
has been used to model CSMA networks [il9i in which nodes are not allowed to transmit if they detect 
an ongoing transmission in their vicinity. 

In the HC-I model we assume that all n interferers in the radius R network are active unless they 
are within a distance h of any other interferer or the representative transmitter. In the HC-II model, 
each interferer is assigned an independent uniform (0, 1] random variable called a mark, where the mark 
assigned to the i-th interferer is Afj. All interferers are active unless they are within a distance h of the 
representative transmitter, or if they are within a distance h of any potential interferer and have a higher 
mark than that potential interferer. 

In comparison, the Matem type I and type II processes are constmcted from an underlying Poisson 
point processes. For the Type I process, each point of the underlying Poisson process that is within a 
distance h from another point of the underlying process is removed. For the Type II process, a point of 
the underlying Poisson process is retained only if there is no other point of the Poisson process within a 
distance h of it which has a mark of lesser value. Thus, these processes are the result of dependent thinning 
of an underlying Poisson point process. To incorporate an active transmitter at a particular location, the 
Matem process has to be conditioned on having a node at that location which leads to added mathematical 
complexity as the probability of a node being at a particular point is zero. We sidestep this complexity 
for HC-I and HC-II by assuming that the representative transmitter is not part of the random process but 
is located at a deterministic point. 

Recall that (Pi, P2, ■ ■ ■ , P„) = g{Xi, X2, ■ ■ ■ , X„; Mi, M2, ■ ■ ■ , Mm', Xt) assigns the transmit powers 
to the interferers and can be defined as follows. For the HC-I model, 

p _ f if |Xj — Xt\ < h or 3j such that \Xi — Xj\ < h 
' ~ \ 1 otherwise. 

For the HC-E model, 

p _ ( if |Xi - Xt\ < h or 3j such that \Xi - Xj\ < h and Mi > Mj 
* 1^ 1 otherwise. 

The HC-I and HC-II models satisfy the asymptotic independence property in (|5), and ([10] ), as shown 
in Appendix |Dl Thus, (fT4l) can be used to approximate the asymptotic and mean bits/symbol for these 
models with the appropriate limiting effective density of transmissions p. The limiting values of p are 
derived in Appendix ID] and are as follows. 
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For the HC-I model, as i? — )■ oo, 

p = pp Pr(Pi = 1) ^ ppexp(-pp7r/i^) . (17) 
For the HC-II model, as i? — )■ oo, 

p^4^(l-exp(-pp7r/i2)). (18) 

V. Frequency Reuse Systems 

In certain wireless communications systems such as cellular networks with Time-Division-Multiple- 
Access (TDMA) and frequency reuse, at most one wireless node is allowed to transmit at a given time, in 
a given frequency band, and in a given spatial region. The active nodes in such a network are of course 
spatially correlated and the framework developed above can be used to analyze such systems. Note that 
cellular systems with regular cells are notoriously difficult to analyze in closed form, even when there 
is universal frequency reuse, as noted in references such as f2(T|. Thus, the asymptotic analysis provided 
here provides insight into a problem that is known to be challenging. 

Consider the network illustrated in Figure [T] which shows a cellular network with reuse factor 3, with 
only one transmitter active in a given time slot, in a given cell. The base stations are represented by the 
dots and the representative transmitter by the square and the in-band transmitters by the crosses. 

To apply the framework developed above, we first assume that the base stations are distributed on 
a plane in a hexagonal grid pattern, with one base station (the representative receiver) located at the 
origin. Let pc denote the area density of the base stations. The representative transmitter is located at a 
deterministic point in the cell at the origin. Overlaid on this grid of base stations is the circular network 
of radius R, centered at the origin with n potential interferers distributed in an i.i.d. fashion, with uniform 
probability in the circular network with their positions denoted as Xi, X2, ■ ■ ■ X„ as in Section III. As in 
the HC-II model, we assign each potential interferers a mark which is an i.i.d. random variable uniformly 
distributed on (0,1]. The mark associated with the i-th potential interferer is denoted by Mj, and is used 
to decide which node transmits in a given cell. Let us assume that the cell at the origin, denoted Co, is 
assigned to frequency band 0, and the other cells assigned to band be denoted Ci, C2, ■ ■ ■ . For simplicity, 
we shall not consider power control in this work (systems which employ power control, and more general 
cell and reuse structures will be treated in a future work). The set of active nodes in band in a given 
time slot is determined by the function g{Xi,X2, ■ ■ ■ , X„; Mi, M2, ■ ■ ■ , Mm] Xt) defined as follows. 

^ _ f 1 if 3 J s.t. Xi G Cj AND Mi < Mk^k^i for which Xk G Cj 
' ~ \ otherwise 

As shown in Appendix [Gl this model satisfies the asymptotic independence property in (|5]), and (fTOl ). 
The limiting effective density p is found in Appendix [G] to be 

p = -pjl-e'-). (20) 



K 

Thus, we can approximate the mean bits/symbol of link with length by (fT4l) and p as above. The 
reuse factor (3 in the example) is denoted by k. Note that if the representative transmitter is located at 
the cell edge, pc = Substituting for pc in (|20l) and the resulting expression for p into ([14)) yields 

the following expression for the mean spectral efficiency of a cell-edge user 



E[C] ^ log. 



1 + 



V 



Ka (It: 
N 7 T-:r sm 

27r^ 1 — e 



a 



(21) 
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Fig. 1. Illustration of a cellular network with frequency reuse of 3. Active mobile nodes in a given time slot at a given frequency band are 
illustrated by the crosses. The representative transmitter is assumed to transmit in the given time-slot and frequency band and is denoted by 
the square, and the base stations are denoted by the dots. The mobile users (both active and inactive) are distributed randomly with uniform 
probability in the circular network of radius R which is then taken to infinity 



When the density of potential interferers to base stations is high, (1211) further simplifies to 



E[C] ^ log, 1 + 



(22) 



Note that the last expression takes a very simple form which depends only on a few parameters. Its 
accuracy over reasonable values of N, a and other parameters is validated in Section I VII I 

VI. Boolean Cluster Networks 

In this section, we illustrate the utility of Theorem 1 to the Boolean cluster model which is described 
in Section I, which also details how it differs from the related Matem cluster model. More precisely, we 
can define the Boolean cluster model as follows. 

Consider a set of cluster centers Mi, M2, ■ ■ • Mm, i.i.d., with uniform probability, in the circular network 
such that m = ixpi^R^, where is the density of the boolean clusters. The function g{Xi, X2, ■ ■ ■ , X„; Mi, M2, 
for this model can be written as follows. 

^_ri ifX,e[jZ^BiM,,h) 

* 1^ otherwise, 

where B(X, d) is a disk of radius d centered on X. Note that this model trivially satisfies the asymptotic 
independence property because as R ^ 00, the probability that two potential interferers Xi and Xj same 
cluster diminishes to zero, making their transmit powers independent asymptotically. The effective density 
of active transmissions is given simply by the product of the density of potential interferers pp and the 
percent coverage of the Boolean model (see e.g. [[T8l ) with circular grains of radius h. This value is 

P = Pp (1 - e-""""') . (24) 

Thus, the mean spectral efficiency is given by (fT4l) with the above definition for the limiting effective 
density of interferers. 

VII. Numerical Simulations 

We conducted Monte Carlo simulations to validate the asymptotic results for the mean bits/symbol 
derived in the previous sections, for the different systems considered in this paper. 
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A. Maternlike hard-core models 

We simulated both HC-I and HC-II networks with a = A. The link-lengths for both sets of simulations 
were selected such that uppr"^ = 1, i.e. on average 1 other potential interferer is closer to the representative 
receiver than the representative transmitter. The radius of the guard zones h were varied as a function of 
rx and are given in the corresponding figures. 

Figure |2] illustrates the mean bits/symbol from 1000 simulations of the HC-I model with guard zones 
that are a multiple of the length of the representative link. Note from the figure that the asymptotic 
approximation is very close to the simulated mean bits/symbol. For > 6 diversity branches, the simulated 
and asymptotic mean bits/symbol (from (fT4l) ') with p defined appropriately) are virtually indistinguishable. 
Figure |3] illustrates the mean bits/symbol from simulations of the HC-II model with guard zones that are 
a multiple of the length of the representative link. Note from the figure that the asymptotic approximation 
is very close to the simulated mean bits/symbol, even when the number of diversity branches is small. 
Figure |4] shows the bits/symbol obtained from 100 simulations of the HC-II process for different receive 
diversity orders. The concentration of points with increasing diversity order indicates that the bits/symbol 
concentrates on the asymptotic prediction which implies convergence in probability. 
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10 20 30 40 50 60 70 

Diversity order N 

Fig. 4. Bits/symbol from 100 simulations of HC-II model of interferers vs. receive diversity order. The link length rr was such that 
■Kp-pV^ = 1, and the radius of the guard zone, h was a multiple of tt as labeled in the figure. 



B. Frequency-reuse Systems 

In order to validate the asymptotic results presented in Section |Vl we simulated cellular networks with 
hexagonal cells and frequency reuse factor 3 as illustrated in Figure [T] The path-loss exponent a = 4, the 
density of potential interferers (wireless nodes) pp = 0.01 and different relative densities of base stations 
to wireless nodes (indicated in the Figure) were used. At most one node was allowed to transmit in any 
active cell. Figure |5] illustrates the mean bits/symbol of the representative link of length r^, where tt 
was such that uppr"^ = 1, i.e. there is on average one potential interferer closer to the representative base 
station than the representative transmitter. 

The markers in Figure |5] illustrate the simulated values and the solid lines represent (fT4] ). with the 
effective density of interferers p given by (|20l) . The asymptotic prediction for the mean spectral efficiencies 
are within 1% of the simulated values for > 6, validating the asymptotic predictions. The dashed lines 
represent the simulated standard deviations which decay with diversity order. For the systems we simulated, 
the standard deviation is less than 5% of the asymptotic predictions when N > 10. The decaying standard 
deviation indicates mean-square convergence which implies convergence in probability. 

Figure |6] represents simulations of a similar network model with Pc/ Pp = 0.05, with the representative 
transmitter located at the edge of the cell centered at the origin. The figure shows the reuse-normalized 
mean bits/symbol vs. the diversity order at the representative base station, i.e. the vertical axis of the graph 
illustrates | times the mean bits per symbol. The simulated values agree with the asymptotic prediction to 
within 1.05 % when the diversity order is six or greater. We also plotted the asymptotic mean bits/second 
for a system with universal frequency reuse on the same graph (dashed line) which illustrates that a several 
fold increase in mean bits/symbol (normalized by reuse factor) is possible by implementing frequency 
reuse, even when the system is interference limited as is assumed throughout this work. 

C. Boolean Cluster Model 

In order to validate the asymptotic results presented in Section |VIl we simulated Boolean cluster 
networks with path-loss exponent a = 4. The density of potential interferers (wireless nodes) pp = 0.01 
and different relative densities of clusters to wireless nodes was used. Figure |7] illustrates the mean 
bits/symbol of the representative link of length tt- tt was such that nppr^ = 1, i.e. there is on average 
one node closer to the representative base station than the representative transmitter. The simulated mean 
bits/symbol are represented by the markers and the asymptotic prediction of (fT4l) with the effective density 
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Fig. 5. Mean bits/symbol vs. receiver diversity order for cellular system with reuse factor of 3 and different relative densities of wireless 
nodes to base stations. Representative link length rr is such that ivppr^ — 1. The three dashed lines (difficult to distinguish on this scale) 
represent the decaying standard deviation of the bits/symbol. 
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Fig. 6. Mean bits/symbol normalized by reuse factor vs. receiver diversity order for cellular system with reuse factor of 3 and relative 
density of base stations to wireless nodes Pc/pp = 0.05. The representative transmitter is located at the edge of the cell. To illustrate, the 
asymptotic mean bits/symbol for a universal frequency reuse system is also plotted. 



of active inteitferers given by (l24T i is represented by the solid line. The simulated mean bits/symbol agrees 
with the asymptotic expression to within 2.5% when the number of diversity branches is 8 or larger in all 
the cases we considered. The dashed lines in the figure indicate the standard deviation of the bits/symbol 
from simulations which decay with increasing diversity order. For the cases we consider, the standard 
deviation is 10% or less of the asymptotic prediction when the number of diversity branches is 12 or 
more. The decay in standard deviation indicates convergence in mean square, and hence in probability. 

VIIL Summary and Conclusions 

An asymptotic technique to compute the mean bits/symbol achievable in wireless networks where active 
interferers are spatially correlated and receivers have N diversity branches with linear MMSE processing 
is introduced. In particular, we consider networks which are the result of the dependent thinning of an 
underlying process with independent node positions. The thinning operation is such that the received 
powers of pairs of nodes satisfy a particular asymptotic independence criteria. This framework is applied 
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Fig. 7. Mean bits/symbol vs. diversity order for Boolean cluster model with varying cluster density. The three dashed lines represent the 
decaying standard deviation of the bits/symbol. 



to two hard-core models, the HC-I and HC-II models defined in Section |IVl which are analogous to 
Matem type I and II models which are frequently used to model networks with CSMA MAC protocols. 
We additionally apply this result to analyze cellular networks with orthogonal transmissions and frequency 
reuse patterns, and a Boolean cluster model which can be used to model networks with hotspots. These 
models, in particular the hard-core and cellular models have proven difficult to analyze in the literature. 

For systems satisfying our model, when the receive diversity order grows large, the influence of the 
spatial dependence of the active nodes on the bits/symbol diminishes. As a result, the mean bits/symbol 
for "large" can be accurately approximated using simple expressions. These results provide insight 
into the benefits of increasing diversity order, which is directly related to cost in many systems, on the 
expected data rate. These results can also be utilized to optimize parameters such as the guard zone- 
radius in hard-core networks, and reuse factors in cellular systems to maximize mean data rates. The 
asymptotic results are validated by Monte carlo simulations which match theoretical predictions even 
when the diversity order is only moderately large. These findings provide insight into wireless networks 
with correlated transmissions, which are well known to be challenging to analyze, and for which only a 
limited number of analytical results exist in the literature. 

Appendix 

A. Proof of Main Result 

The main property used in proving Theorem 1 is the following lemma which is a variant of Lemma 1 
from Govindasamy et al. [13] and Lemma 4.3 of Tse and Hanly [|2TI . 
Lemma 1: Consider the quantity 

^" = ^^'(^^8*8*)'^ (25) 

where s G C^^^ and S G C^^" comprise i.i.d., zero-mean, unit-variance entries, n/N = c, and ^' = 
diag(?/^i, ^p2, ■ ■ ■ ipn)- Suppose that as n, N ^ oo, the e.d.f. of (^^i, ^"2, ■ ■ ■ i^n) converges in probability to 
H{x). Additionally, assume that there exists an such that Vn > hq ,Xmin > ^^ib > 0, with 

probability 1, where Amj„(A) denotes the minimum eigenvalue of A. Then, 77V — 7 in probability where 
7 equals the non-negative real solution for /3 in ([8]). 

Proof: By assumption, the e.d.f. of the diagonal entries of converges to a non-random limit in 
probability. Hence, by Theorem 4.1 of Bai and Silverstein [|22|. the e.d.f of the eigenvalues of -^S^'S^ 



13 



V 

n 



> f ) -> 0, which 



converges in probability to a non-random limit. From Lemma 4.3 of [1211 and Lemma 1 of [fT3l . it is 
evident that Lemma 1 holds if the e.d.f. of the eigenvalues of -^S^'S^ converges in probability to a 
non-random limit as n, — )■ oo. ■ 
Up to a scale factor, (l25l) has the same form as the equation for the normalized SIR /^at in dS]) with 
s = gr, * = diag [m/^Pir^'^, A^^/^Pgr^", ■ ■ -, m/^P^r-'^), and S = [gi g2 ■ ■ ■ g„]. To use this resuk 
to characterize (3^, v^e first condition on the event (denoted by V) that the fraction of potential interferers 
that are active is greater than | of its limiting value. This conditioning assures that the minimum eigenvalue 
condition requierd in the lemma is satisfied. Note that the number | here is somewhat arbitrary as any 
positive rational number less than one would suffice. Thus, V is the event that ^ > ,. Since ^ — )■ z/ 

in probability by the asymptotic independence property, as n, N, R ^ oo,Pt ^ — 
implies that 

Pr(I?) = Pr fffl > — ^ ^ 1. (26) 
\ n A J 

We shall use this fact to remove the conditioning on V. 

To use Lemma [H we first show that conditioned on V, the e.d.f. of the pj„ = N°'^'^Pir^" terms, defined 
in (|7]), converges in probability to a limiting function H(x) in Lemma [2l We then show that the minimum 
eigenvalue property is satisfied in Lemma |3] The proof is completed by evaluating ([8]) with H(x) from 
(1281) in Lemma H 

Lemma 2: If dH) holds and conditioned on V, as n, N, R oo as in Theorem 1, 

HJx) Hix) = lim - V Pr(P, = 1) Pr (A^tr"" < xlPi = l) + - V Pr(Pi = 0) , (27) 

i=l 1=1 



where the convergence is in probability. If (|TOl ) holds, then 

Proof: Given in Appendix 

Lemma 3: Let ^ and S be as defined above. Conditioned on V there exists an no such that for all 
n > uq, and with probability 1, 

A„..„ (ls*St) . > (^)' ^ (l - . (29, 

Proof: Given in Appendix \C\ 

Thus, from Lemmas \T\ |2] and |3l conditioned on V, (3^ ^ /3 in probability, i.e. for any e > 

Pii\f3N-f3\>e\V)^0. (30) 
We remove the conditioning by using (|30] ) and (|26] ) as follows. For any e > 0, 

Pr - /3| > e) = Pr {{P^ - /3| > e |I)) Pr {V) + Pr (1/3^ - f3\ > e |P^) Pr {V') ^ 0. (31) 

The first part of Theorem 1 is thus proved. If ([TO] ) holds, we can substitute the derivative of H{x) from 
(l28l) . into ([8]), and integrate to get 

CSC — - ——^2Fi 1,1 ,2 ; = 1 (32) 



\ J {a — 2)c^~Q \ a a 

This integral (without the scale factor of lim„^oo Pr(-Pi = 1)) has been evaluated before in [17.1 , and we 
refer the reader there for the details of its computation. Applying the Pfaff transform (see e.g. [fT6l ) to 
the hypergeometric function in (|32] ) yields (fTTT) . 
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B. Proof of Lemma H] 

We first will show that this property holds without the conditioning. Recalling the definition of ain{x) 
from Section Hill note that E[ain{x)\ = Pr(pm < x) and E[ain{x)ajn{x)\ = Pr(pj„ < x.pjn < x). Next, 
define the following 



^ n 1 

Hnix) = - V^K(x)] = - VPr(rr"P,iV-/2 < x) 



H^(x) = lim HJx) 



(33) 
(34) 



Recall the definition of Hn{x) from Section Hill Since (|34l) converges, for each e > 0, 3na such that 

for all n > Ua, \Hn{x) — HaQ{x)\ < |e. Thus Vn > Ua, 



Ft {\H^{x) ~ H^{x)\ >e) < Pr - | >| 



Pr 



i=l 
n n 



i=l 



ne , 



var <^ ^ai„(x) 



i=l 



J^^cov(ai„(x),aj„(x)) = ^-^ ^ ^(E[ai„(x)aj„(x)] - E[ai„(x)]E[aj„(x)]) 

i=l j=l i=l j=l 

n n 

^ ^ (Pr(pi„ < x,pj>, < x) - Pr(pi„ < x)Pr(pj„ < x)) 
i=i j=i 



(35) 



where the inequality follows from the Chebyschev inequality. As n, N, R ^ oo, from (HJ we have (|35T ) 
going to zero. Now, if we let 



H{x) = Hoo{x) , 



(36) 



Pr {\Hn{x) — H{x)\ > e) — 7- which implies that Hn{x) converges in probability to H{x). Now since V 
holds with probability approaching unity, we have 

Pr {\Hnix) - H{x)\ >e\V) = [Pr (|i/„(x) - H{x)\ > e) - 



Pr (V) 

Pr {\Hn{x) - H{x)\ > e \V')Ft (V)] 



(37) 



To show (1271) . we note from (1331) . (1361) . that as n,N,R ^ oo. 



H(x) 



1 " 

- ^[Pr {PiN'^r-'^ < x|P, = l)Pr(Pi = 1) + Pr [P.N^rr'^ < x\Pi = 0)Pr(P, = 0)] 

^ i=l 

which implies (l27l) since Pr {^PiN^r~°' < x\Pi = O) = 1. 

For each n, the rj's are identically distributed. Thus, if (ITOl ) holds, (1271 ) becomes 

if(a;) = lim - V Pr(Pi = 1) Pr (A^trf" < x) + lim - V Pr(Pi = 0) 



Additionally, we note from lIH that the CDF of A^"/Vf is 

Pr(Ar"/Vr" < x) = 1 - ^x'^l 



i=l 



(38) 



(39) 



15 



which does not depend on N because A^, n, and R are related such that the CDF of N°'/'^r^ ° is independent 
of A^. Substituting dlH) into dMI) yields the following as N,n,R^ oo, 



H(x) 



lim -VPr(P, = 1) I (l-^x-^l,^,,../2^ ^ + lim - VPr(P, = 0). 

\ i=l / ^ ' i=l 



(1281) follows from substituting p = (lim„^oo ^ Ya=i ^^i^i = 1)) Pp and the fact that Pr(Pi = l) + Pr(Pi 
0) = 1, into the previous expression. 

C. Proof of Lemma \3\ 

Recalling that T is the set of active transmitters, write (|29] ) as 



n n 
i=l i=l 

= E g'si + - R-ns^s^ . (40) 

Let lim„_>.oo ^ ~^ which holds in probability from the asymptotic independence property described in 
Section II. u represents the fraction of potential interferers that are active in the limit. Let the members of 
T be labeled ii,i2, - ■ ■ . Next, define a new set of indices T' which contains the indices of at most ^un 
of the active nodes (note that any positive fraction less than one, would work here). In other words, 

r = lik:k< mini ^iyn,\T\]\ (41) 



.4 

Since c = n/N and n = irppR^, R = \J Thus, (l40l) can be written as 

' ' ieT' ier\T' ieT 

For notational simplicity, define the following matrices 

'TCPp\ f c\T'\ 1 



ieT' 



/Yf-i 

= E + E("«"" - ^~")&gl • (43) 

ieT\T' ieT 

Since T2 is non-negative definite, Am,m(T2) > 0. Thus, by the Weyl inequality (see e.g. [|23l ). 



A. 



mm(Ti) -|- Amm(T2) ^ ^min (Ti) . (44) 



Conditioned on "D, ^ = ^, and from Theorem 5.11 of [22], with probability 1, 



Substituting into (142)) . 

/ Q / / \ ^ 

/ /7r/)„\ # .Sn/y / / 4 \ 

P = 1 (46) 



Pr I lim X^iJT 
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The proof is completed by combining (|46] ). P4l) and Theorem 6.3 of [|22]| . which yields that conditioned 
on there exists an uq for which for all n > uq, with probability 1, 

A.„.„ (1S*S') > A„„(TO > (^) ' !^ j^l - \ (47) 

D. Proof that HC-I and HC-II models satisfy the asymptotic independence property 




Representative 
transmitter 



Fig. 8. Hard-core network illustrating the representative receiver, representative transmitter and two interferers. 

In this section, we show that both the HC-I and HC-II models satisfy ([5]). Instead of showing that 
the joint CDFs factor in the limit, it will be simpler to show that the joint Complementary CDFs (i.e. 
1-CDF) factor instead. Figure [8] illustrates the network with the representative receiver at the origin, the 
representative transmitter at a distance from the origin, and interferers i and j. Let A denote the shaded 
region which is an annulus centered at the origin with inner radius tt + 2h, and outer radius R — h. Let 
A represent the event that Xj G A,Xj G A and |Xj — Xj\ > 2h. We consider Xi,Xj G A in order to 
eliminate edge effects and the effects of the representative transmitter. Note that as R ^ oo, Pr(^) — t- 1 
since Pr(Xj G A, Xj G A) — )• 1 as the fraction of the area of the circle of radius R that is not in A goes 
to zero, and the probability that the distance between two random points in A is less than 2h becomes 
negligible since this latter probability is of order 0(-^) (e-g- see [|24l Equation 2.3.63). 

Since Fr{PiN"/^r;'' > x, P^N^'/^rJ'' > x) > only if P, = Pj = 1, we have the following 

Pr(PiA^tr-" > x.PjN'^r-'^ > x) 

= Pr(iVtr-° > X, N^r-'^ > x\A, Pi = 1, Pj = 1) Pr(Pi = 1, Pj = 1\A) Pr(^) 

+ Pr(A^trr" > x, mrf^ > xlA", Pi = 1, Pj = 1) Pr(Pi = 1, Pj = llA") Ft{A'). (48) 

Since as i? — > oo, Pr(^'^) 0, we only need to consider the first term on the RHS of (|48T) . The limit of 
this term can be found using the following two lemmas. 

Lemma 4: For both the HC-I and HC-II models, and i ^ j, we have 

lim Pr(A^tr-" > x, A^trT" > x\A, Pi = 1,P = 1) 

= lim Pr(A^^r"" > x) lim FT{N^r~'^ > x) and, (49) 
lim Pr(A^trr" > x, \A, Pi = 1) = lim Pr(A^tr,-" > x) (50) 

n— >oo n— >oo 
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Proof: Given in Appendix E 

Lemma 5: For both the HC-I and HC-II models, and i ^ j, 

lim Pr(Pi = l,Pj = 1\A) = lim Pr(Pi = 1) lim Pr(Pj = 1) (51) 



n—^oo 
1.2 



where lim„^oo Pr(Pi = 1) = e'^P^^^ and lim^^oo Pr(Pi = 1) = for the HC-I and HC-II models 

respectively. 

Proof: Given in Appendix F 

Substituting Pr(j4) 1, dSB and (|49l) into ([48]), and taking the limit, we have for i ^ j, 

lim PY^PiN^r-'' > x,PjN'^rj" > x) 

= lim Pr(Pj = 1) lim Pr(P,- = 1) lim Pr(iV^r"" > x) lim Pr(A^^r~" > x) 

= lim Pr(P,iV^r-" > x) lim Pr(PiA^trr" > x) , (52) 

n— >oo n— ^oo 

where §^ is from (l49]) and the fact that Pr(Pj = 1|^) ^ Pr(Pj = 1). Recalling that pin = PiN^r~°' 
and pjn = PjN'^r-"', ^ implies ©. Additionally, Pr(^) 1 and ^ imply (fTOl) . 

£. Proof of Lemma E] 

Recall that ^ is the event that \Xi -Xj\ > 2h and Xi, Xj e A and let \A\ = tt{R - - n{rT + 2h)^ 
be the area of the region A. Thus, 

Pr(A^f r"" > X, N'^rj'' > x\A, Pi = 1, Pj = 1) 

= Pr(A^trr° > x\A, N'^rJ'' > x) Pr(Artr-" > x\A) (53) 
Noting that A implies that X^ must be outside B{Xi, 2h), 

Pr(A^trr° > x\A,N'^rJ°' > x) = Pr(ri < y/Nx'^\A,rj < \fNx~^) 
ttNx~^ \ 



< 



vr (^p2) ^-^ 



vr(p - /,)2 _ Jr(rT + 2/i)2 - 47r/i2 j ^{JVT(ij-h)-<x<ivf (r^+2h)-} + l{x<ivf 
where we have used n/N = c and ©. Taking the limit as n, N, R ^ oo and applying ([T]) again. 



2 



Similarly, we can write a lower bound which is active if B{Xj, 2h) C P(0, yNx «) as follows 

Pr(A^trr" > a;|^,iVtrj" > a;) = Pr(ri < ^x~^\A,rj < ViVx"^) 



> 



_ 47r/;,2 J ^ {n^ (R-h)-°'<x<N^ {rT+2h)-°'} + {a;<Af f } 

Comparing with (|54l ) and recalling that 1^41 = 0(P^) = 0(A^), note that the lower bound converges to 
the same form as the upper bound as n, N, R ^ oo. Thus, 

Since Xi and Xj are independent and Pr(^) — )■ 1, from (|55T ). 

lim Pr(A^tr7° > xU) = lim PrfA^^r"" > x) = (—) x'^lr,^^ ^ + 1 

71— >00 -' 71— ^OO V C / {{-^] <x\ < 



Taking the limit of (|53l) , substituting (1551 ) and the previous expression yields (|49l ). (|50l) follows directly 
since conditioned on ^, Pj is independent of the specific location of Xi, and hence rj. 
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F. Proof of Lemma \5\ 

Let #B{Xi,h) = \{i : i ^ i, Xi E B{Xi,h)}\ denote the number of potential interferers in B{Xi,h) 
excluding the interferer at Xj itself. For both the HC-I and HC-II models, conditioned on #B{Xi, h) = k 
and #B{Xj, h) = m, the events Pi = 1 and Pj = 1 are independent. Thus, 

Pr(P, = l,Pj = l,#B{Xi,h) = k,#B{Xj,h) = m\A) = PT{#B{Xi,h) = k,#B{Xj,h) = m\A) 

X Pr(Pi = l\A,#B{X„h) = k)Fi{Pj = l\A,#B{Xj,h) = m). (56) 

Summing (|56l) over all m and k excluding the potential interferers at Xi and Xj yields, 

n-2 n-2 

Pr(Pi = l,Pj = 1\A) = ^^Pr(Pi = l\A,#B{X,,h) = k) Pr(P,- = l\A,#B{X^,h) = m) 

k=0 m=0 

Pr(#P(X„ h) = k,#B{Xj, h) = m\A). (57) 

For the HC-I model. Pi = I and Pj = 1 if and only if #B(Xi, h) = and #B(Xj , h) = 0, respectively. 
Hence (1571) becomes 

Pr(P, = l,Pj = 1|^) = FT{#B{Xi,h) = 0,#B{Xj,h) = 0\A). 

= Pr(#P(Xi, h) = 0\A,#B{Xj, h) = 0) FT{#B{Xj, h) = 0,^) . (58) 

Recall that conditioned on A, Xi ^ B(Xj, h), and Pj = 1 if the n — 2 remaining potential interferers are 
not in B{Xj, h). Using the fact that n = 'KppE? we have 

lim Pr(#P(X,-, h) = 0\A) = lim 1 — = e"""'''^ . (59) 

n-s>oo R^oo y TtR^ J 

Furthermore, since Pr(^) — i- 1 as n, X, P — > oo, we have 

lim Pr(P,- = 1) = lim [Fr{#B{X^, h) = 0\A) Fr{A) + Pr(P,- = 1, A^)] = e"^''^^' . (60) 

Conditioned on #B{Xj, h) = the n — 2 potential interferers are outside B(Xj, h). Since Pj = 1 only if 
the remaining n — 2 potential interferers are outside B(Xi, h) as well, 

Pr(#P(X„ h) = 0\A,#B{Xj, h)=0)= [1- 



TTR^-TTh\ 

lim Pr(#P(Xi, h) = 0\A,#B{Xj, h) = 0) = e'^P^''^ = lim Pr(Pi = 1) (61) 



where the last equality follows from symmetry between Xj and Xj and (|60l) . Substituting (1601) and (|6T 
into (I58l ) and taking the limit yields dlB for the HC-I model. 

For the HC-II model, using the fact that Pr(Pj = l|#P(Xi, h) = k,A) = dSH) becomes 



n-2 n-2 



Pr(P, = 1,P,- = 1|^) = ^ ^ ^^^Pr(#5(X„/.) = k,#B{X„h) = m\A), 

n—2 ^ n—k—2 _ 

V Fii#BiXi,h) = k\A) V Pr(#P(X,-,/i) =m|#P(Xj,/i) = fc,^), 

^-^ k + 1 m + 1 



fe=0 ' m=0 



^ A;+ 1 V k \~^) 
^ m + 1 V ni ) \7fR^ -TThy V 7rP2 _ T,h^ ) ^ 



k=0 

n—k—2 _ / , „\ / , o \ m / 79 \ n— fc— 2— m 
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Consider the previous expression written in terms of a double series (please see e.g. [|25l for basic 
properties of double series) with the variables rii and n2 replacing n, in the terms in the summation, and 
the upper limit of the summation respectively. Additionally, since R = 



n 
wpp ' 



n2-2 



k=0 



h 



2 \ 



El ( rii — k — 2 
m + 1 



m=0 



m 



+ l\ k J \ni/pp 



2 \ ni~k—2 



X 



ni/pp - 



1 



rii/pp - 



ni — k—2—m 



Taking the limits of (|63l) as ni — oo followed by n2 oo, we have 



lim lim Pr(Pi = 1, R = 1|^) = lim V 



n2— ^-oo ni— >oo 



lim lim Pr(Pi = = 1|^) 

n2— >oo ni— >oo 



n2-s>oo ^ — ' {k 
k=0 ^ 



1)! S 



(7r/)6/i2)"^e-'^'''''f 



m=0 



(m + 1)! 



Since both the sums on the RHS of (1641) are finite for all n2, by Theorem 2.11 of [|25l . 



(63) 

(64) 
(65) 



lim Pr(Pi = 1, P. = 1|^) = lim Pr(Pi = 1, P,- = 1|^) 

n— >oo ni,n2— >c« 



= lim lim Pr(Pi = 1, P- = 1|^) = 
Following similar steps, we have 

lim Pr(Pi = 1) = lim Pr(Pj = 1) 

n—^oo n—^oo 

Substituting (l67l) into (l65l) yields (|5TI) completing the proof. 



npbh'^ 



(66) 



(67) 



G. Proof that Frequency Reuse Model Satisfies Asymptotic Independence Criteria 

Let i3 denote the event that Xi and Xj are in cells assigned to band which are wholly contained in 
B{0,R), and Xi and Xj are in different cells. The event B can be described as: 3k,£ such that Xi E 
Ck,Xj e Ci,Ck C B{0,R),Ce C P(0,P) and k ^ i. Note that Pr(i3) as P ^ oo, since the 

probability that Xi or Xj are in a cell assigned to band zero approaches and the probability that they 
are in the same cell goes to zero. Let Ca denote the union of all cells wholly contained in P(0,P), i.e. 
B{Q,R) ^k- Now consider the following 



Ca n P(0, \ 






Col 




Ca 




Co 





< Pr ( p. "Ti < VNx-^\Pj "Vj < VNx-^,B, Pi = 1, Pj = 1 



< 



Ca n P(0, V 






Ca 




Co 





(68) 



The lower bound holds when the cell containing Xj is a subset of P(0, v Nx"^). The upper bound holds 
if the cell containing Xj intersects with P(0, \fNx~'^) at a single point Xj. As P — )■ oo, the edge effects 
diminish, and the upper and lower bounds converge to same form: 



lim 

n— >oo 



Ca n P(0, ^Nx 



\c. 



— X a 1 



Pr(r,- < \fNx a ' 



(69) 
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which implies that 

lim Pr (P^"ri < VNx'^\B, P^"r,j < ^Nx'^^Pi = 1, Pj = l) = Pr (vi < ViVx"- V (70) 

Following similar steps, we can show the following 

lim Pr (P^"n < VNx-^\B, Pi = 1, Pj = l) = Pr(rj- < VNx'^) (71) 

lim Pr (P,~^ri < y/Nx'^lXi G Ca, Pi = l] = Pr(ri < VNx'^) (72) 

Since as n,N,R ^ oo, Pr(^|Pi = 1) ^ 1, (|72l) implies dTO]). Recalling that if B holds, Xi and are 
in different cells and the probability that Pj = 1 if there are k additional potential interferers in the cell 
containing Xi is we have the following 

Pr(Pi = l,Pj = l\B,m additional potential interferers in cell containing Xj) 

~m + l^A; + lV k )\nR^-\Co\) \ 7rR^-\Co\) 
Multiplying by Pr(mpotential interferers in the cell containing Xj\B), and summing over m, 

n—2 / „\ / 1^ I \ m / 1^ I \ n—2—m 



m=0 

2 k j[nR^-\Co\) V~ 7rR^-\Co\) ^^^^ 



Note the similarity of (1731 ) and (|62l) . Following steps used to derive (|66l) . it can be shown that 

1 _ e-|Co|pp^ 2 



lim Pr(P, = 1, P,- = 1\B)=[ ——. (74) 

]^ _ p-\Co\pp 

lim Pr(P, = 1, \X, G Ca) = (75) 

Taking the product of (1721) and (fTST ). substituting (|69] ), and using the fact that Pr(Xj e Ca) — ^ ^, 

lim Pr {P,r-"m < x) = - [ i^x-A ^ (76) 



As P — )■ OO, Pr(Pj = 1, Pj = 1\B^) — )■ 0, because the probability that Xi and Xj are in the same cell 
— )• 1, which means that in the limit, B'^ holds only when Xi or Xj are not in cells assigned to band 0. 
Taking the product of dTO]), (EB, dH]) and lim^^oo Pr(i3) = yields, 

1 /vro 2\2 / 1 — e^l''°l''p \ ^ 
lim Pr (Pir-°X2 > x, Pr^^X^ > a;) = -^x"- — — (77) 

n->oo ^ V C / \ \Cq\Pp J 

= lim Pr iPirr'^N'^ > x) lim Pr (Pr^^Xt > x) . (78) 

Recalling the definitions of Pin and (|78] ) implies ([5]). Substituting the density of cells pc = into 
(fTSl) and observing that Pr(Xj G (^a) — ^ yields the effective density p in (|20l) . 
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